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We study the generation of terahertz radiation from atoms and molecules driven by an ultrashort 
fundamental laser and its second harmonic field by solving time-dependent Schrodinger equation 
(TDSE). The comparisons between one-, two-, and three- dimensional TDSE numerical simula¬ 
tions show that initial ionized wave-packet and its subsequent acceleration in the laser field and 
rescattering with long-range Coulomb potential play key roles. We also present the dependence of 
the optimum phase delay and yield of terahertz radiation on the laser intensity, wavelength, dura¬ 
tion, and the ratio of two-color laser components. Terahertz wave generation from model hydrogen 
molecules are further investigated by comparing with high harmonic emission. It is found that the 
terahertz yield is following the alignment dependence of ionization rate, while the optimal two-color 
phase delays varies by a small amount when the alignment angle changes from 0 to 90 degrees, 
which reflects alignment dependence of attosecond electron dynamics. Finally we show that tera¬ 
hertz emission might be used to clarify the origin of interference in high harmonic generation from 
aligned molecules by coincidently measuring the angle-resolved THz yields. 

PACS numbers: 33.20.Xx, 42.50.Gy, 42.65.Ky 


I. INTRODUCTION 

Terahertz (THz) radiation, with wavelength between 
0.03 to 3 millimeters, has been found useful for appli¬ 
cations in information and communication technology, 
homeland security, global environmental monitoring and 
ultrafast computing, to name a few [1, 2]. Scientifically, 
it provides resonant access into and hence enables the 
probing of various of modes, such as the motion of free 
electrons, the rotation of molecules, the vibration of crys¬ 
tal lattices and the precessing of spins [3]. With power 
continually being increased and duration being shorten 
into single or half cycle, pulsed intense THz sources are 
now capable of driving and steering non-resonant ultra¬ 
fast processes in matter in unique ways. New phenom¬ 
ena are being observed, e.g., THz harmonic generation 
[4], strong-field control of spin excitation [5] and phase- 
transition [6], etc. 

The efficient generation of THz wave, however, is found 
still challenging. The popular approaches include op¬ 
tical THz generation, solid-state electronics and quan¬ 
tum cascade laser, which are limited either by low fre¬ 
quency conversion efficiency, low temperature require¬ 
ment, or the lack of eigen oscillator [7]. Among these 
techniques, laser air photonics is capable of generating 
THz field strength greater than 1 MV/cm with band¬ 
width of over 100 THz through plasma formation by laser 
ionized gaseous medium [8] . But the mechanism of THz 
generation from the laser induced plasma turns out be 
complicated [9] involved of wave propagation, plasma 
formation, oscillation and collision. For example, the 
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THz can be emitted by the radical acceleration of the 
ionized electrons due to the ponderomotive force gener¬ 
ated by the radial intensity gradient of the optical beam 
[10], or by linear mode conversion from the laser-induced 
plasma wakefield in inhomogeneous plasma [11, 12], or by 
Cerenkov-type mechanism from light propagating within 
the filament [13]. 

In this work, we restrict our study to single atom re¬ 
sponse by focusing on THz generation from gases of much 
lower densities in low pressure and minimizing the collec¬ 
tive motion of plasma and its modification on the light 
propagation. Even in this much simpler case, the mech¬ 
anism of THz emission is still under debated, partly be¬ 
cause of the involved complex dynamics during laser in¬ 
teraction with gases. Nevertheless, down to the funda¬ 
mental origin of THz emission, it could only either be 
the induced current or the induced polarization. There¬ 
fore two popular mechanism are proposed, a microscopic 
photocurrent (PC) model with tunneling photoionization 
[14, 15] and a third-order nonlinear model with four-wave 
mixing (FWM) [16-18]. 

Both models qualitatively agree with the experimental 
findings that the broken symmetry of the laser-gas in¬ 
teraction could enhance the field strength of the emitted 
THz waves by several order-of-magnitudes, when either 
applying a DC bias [19, 20], an AC-bias by using two- 
color pulses [16] or using carrier-envelop phase stabilized 
few-cycle pulse [21]. As a representative case, THz emis¬ 
sion in two-color laser fields has been investigated more 
intensively [14-17, 21-28]. So far, more evidences sup¬ 
port the model based on photocurrent as it is found ion¬ 
ization is indispensable [16], although the physics behind 
photocurrent formation needs to be explored carefully 
[29]. However, in view of the success of four-wave mix¬ 
ing model in applying to the THz detection through air- 
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breakdown coherent detection scheme [18, 30], the two 
models must be closely connected. 

The full quantum description in our previous work and 
others [29, 31-33] provides more insight into the mech¬ 
anism of THz emission. By connecting THz emission 
with high harmonic generation (HHG) in two-color laser 
pulses, the dynamics responsible for THz emission is at¬ 
tributed to the continuum-continuum transition during 
laser-assisted soft collision of the ionized electron with 
the atomic core, while the hard collision with the atomic 
core leads to HHG upon recombination [29, 30]. 

Because the fundamental pulse is strong enough to 
induce ionization, its interaction with the gas can not 
be treated perturbatively. The interaction with the sec¬ 
ond harmonic pulse is however weak which seems to be 
possibly treated perturbatively on the top of the non- 
perturbative behavior from the fundamental pulse. In 
this sense, the photocurrent model can be understood 
by another type of four-wave mixing model, but differed 
from the usual case that the involved lower and up states 
could be really populated by multiphoton ionization or 
tunneling ionization. Therefore the two models can be 
unified as hinted in [9, 29-31]. 

Although the second harmonic pulse is very weak, it 
serves as a gate controlling the phase of the electron 
wavepacket and hence breaks the left to right symmetry 
leading to the formation of net current and the emission 
of THz waves. At the same time, the out of phase high 
harmonic emission between adjacent half cycle in single 
color pulse is changed as well leading to the emission of 
even-order high harmonics [34-38]. Therefore, by con¬ 
trolling the phase delay between the fundamental and 
the second harmonic pulse, the THz and high harmonic 
yield can be controlled simultaneously. This enables the 
calibration of the optimal phase delay (OPD) maximum 
THz yields and it is found that the OPDs related to the 
atomic potential [29, 30]. It is thus expected that the 
sensitivity of THz yields to the Coulomb potential might 
be used to map atomic fields from within and to help the 
full characterization of the rescattering wave packet [29] . 

In the present investigation, we explore further the ef¬ 
fects of atomic potential and other laser parameters on 
the generation of THz for both a hydrogen atom and 
a model Hj molecule interacting with two-color pulses. 
The paper is organized as follows. In section II, we 
present the theory and model for THz generation from 
ionizing atoms by two-color laser pulses. In section III, 
we show the calculated results and discuss the physics 
behind. Finally we draw conclusions in section IV. 


II. THEORY AND METHODS 
A. Photocurrent model 

Under strong laser fields, the electron current is pre¬ 
sumably contributed by photoelectrons released after ion¬ 
ization. When the laser photon energy is much less than 


the ionization potential I p , and I p < 2U P where the pon- 
deromotive energy U p is the averaged kinetic energy of an 
electron moving in the laser field, the mechanism of ion¬ 
ization can be considered as tunneling of electron through 
the potential barrier formed by the atomic potential and 
the instantaneous electric dipole potential. The instant 
ionization rate in the electric field of E(t) can be obtained 
using Ammosov-Delone-Krainov (ADK) model [39, 40] 

!!>(()= Ip |C„. | 3 «P{—5^} (1) 

where Z c and I p are the net resulting charge and ion¬ 
ization potential of the atom, I p h = 0.5au is the ioniza¬ 
tion potential of the hydrogen atom, respectively. And 

k = n* = Z^Iph/Ip, C n * = „. r( J +1)r( „. ) - 
If the laser field strength is much higher such that the 
binding energy of the electron is higher than the height of 
the potential barrier, ionization can be proceed through 
over-the-barrier ionization (OTBI). In OTBI range, we 
can use the following empirical formula suggested by [41], 

WQTBiit) = w TI (t ) (2) 

Once the instantaneous ionization rate is found, the 
free electron density as a function of time can be obtained 
by 


; (t) = n 0 < 1 - exp 


)(t') dt' 


( 3 ) 


where n o is the neutral atom density. Assuming the sub¬ 
sequent motion of the electron after ionization is domi¬ 
nated by the laser field, the photon-current and its accel¬ 
eration contributed by the ionization event U are 


j(U,t) = - n 0 (ti)w(ti) 


p(U) + E(t') dt' 


a(ti,t ) = - n 0 (ti) w(U) E(t). 


Note that we take the initial longitudinal momentum 
p(ti ) = 0 at the instant ionization. The total photon- 
current J(f) and dipole acceleration A(t) are thus formu¬ 
lated as 


j(ti,t)dti, 

J — OO 

A(t) = f a(ti, t) dti. 

J —OO 

From the current, the net residual current density 
(RCD) after the pulse is given by Jrcd = J(t —> oo). 
The frequency-domain spectrum F( fl) are calculated by 
the Fourier transformation of acceleration A(t) with in¬ 
tensity given by G(fl) = |F(fl)| 2 . The total yields of 
THz radiation is given by the integral of the spectra in¬ 
tensity from 0 up to 30 THz. 
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B. Time-dependent Schrodinger equation 


The time-dependent Schrodinger equation (TDSE) for 
an atom interaction with the laser field in dipole approx¬ 
imation is given by 


I"a2 


i s * ( r,f) = 


+ V (r) + r • E(t) 


T(r,f) (4) 


where the linearly polarized (along z direction) two-color 
laser field is given by 



—Ionization 
lOTHz 
30THz 
RCD 
RCDPC 


E(t) = [E u cos (wt) + E 2u] cos (2wt + cj>)] zf{t). (5) 

Here, E u , E 2u , f(t), u>, 4> denote fundamental and 
its second-harmonic laser field strength, temporal en¬ 
velop, fundamental frequency, two-color phase delay, re¬ 
spectively. By numerically solving the TDSE in three- 
dimension, the induced dipole acceleration along z direc¬ 
tion is given by the quantum average: 

a(t) = (^-E(t)), (6) 

from which the THz radiation can be obtained by Fourier 
transformation. The corresponding electron current den¬ 
sity [42] is given by J(f) = ( t ) dt and can be recast 

into 

J(t) = J'(t) + J bb (t) 

where the fast oscillation part Jbbit) results from the 
dipole transitions between the populated Rydberg states, 
and J' it) is the rest part of the current, which generates 
the broadband terahertz radiation. The residual current 
density can be found numerically by time-averaging the 
current after the laser turns off to cancel the Jbbit) com¬ 
ponent, Jrgd = 7 ftl +T dt J{t) which can be directly 
compared to that from the PC model. The THz spectra 
and yields can be computed from the dipole acceleration 
at the same manner given in the previous subsection. 

As an example, we show in Fig. 1 the typical THz yield 
modulations from an Hydrogen atom in the two-color 
laser pulse, obtained from the TDSE simulation. As the 
two-color phase delay varies, the modulations of terahertz 
yield below lOTHz, 30THz and RCD takes the maximum 
at nearly the same optimum phase delay (OPD) (j) m as 
shown in Fig. 1. It suggests that the RCD can be used as 
a good parameter characterizing the phase-delay depen¬ 
dence of THz yields. On the other hand, the ionization 
probability is less modulated by the applied weak second 
harmonic pulse which indicates that the modulation of 
THz yield is mainly contributed by the controlling of the 
electron dynamics after ionization in accordance to the 
analysis presented in [29, 30]. For comparison, we show 
the modulation of RCD from PC model, which deviates 
from the TDSE calculation because of the neglecting of 
the Coulomb potential in the dynamics of free electrons 
(???). 


FIG. 1. (color online) Modulation of THz yields, residual 
current density and ionization probability as functions of the 
phase-delay between the two color pulses. The TDSE calcula¬ 
tion is carried out for hydrogen atoms. The fundamental laser 
pulse has the intensity of 2 x 10 14 W/cm 2 , the wavelength of 
800nm, the FWHM of 25fs, and the second harmonic intensity 
is 1% of fundamental pulse. 

III. RESULTS AND DISCUSSIONS 
A. Laser intensity, wavelength and atomic orbitals 

We first examine in details how the THz yields from 
hydrogen atoms depend on various of parameters of the 
two-color laser pulses, especially we pay attention to the 
optimal phase-delay that maximizes the THz yields. In 
Fig. 2, the laser-intensity dependence of the optimal 
phases and the corresponding THz yields are compared 
between different models. The ionization potential is 
adjusted taking the same value of I p = 0.5au. When 
the atomic potential is the long-range Coulomb poten¬ 
tial, for either ID, 2D or 3D calculations, the optimal 
phase follows the same trend that changes from 0 . 97 T to 
0 . 67 T as the laser intensity varies from Iq to 10/o where 
Io = 10 14 W/cm 2 . In contrast, the optimal phase pre¬ 
dicted by the PC model varies little around 0 . 57 T (the 
variation might be due to the depletion of the ground 
state). It confirms our previous investigation that PC 
model for THz generation fails due to the ignorance of 
the Coulomb potential [29]. When a short range poten¬ 
tial rather than the Coulomb potential is used, such as 
for the negative ions, the optimal phase of THz genera¬ 
tion shows complete different behavior, increasing from 
0.17T to 0 . 57 T as the laser intensity increases as shown in 
Fig. 2(a). It clearly demonstrates the importance of the 
long range potential. 

The variation of the OPD with laser intensity indicates 
that different ionization mechanism are involved. Let us 
now focus on the optimal phases from the 3D calcula¬ 
tions with Coulomb potentials that are closely related 
to the experimental observation. As the laser intensity 
increases, the ionization mechanisms changes from multi- 
photon ionization (MPI) to tunneling ionization (TI) and 
then over-the-barrier ionization (OTBI). 

In the regime of MPI when 7 > 1, terahertz genera¬ 
tion is dominated by four-wave-mixing [31] with predic- 
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FIG. 2. (color online) The optimal yields (a) and the corre¬ 
sponding optimal phase-delay (b) of THz generation at vari¬ 
ous of laser intensities calculated by ID- (black circle), by 2D- 
(red square), by 3D- (green diamond) TDSE with Coulomb 
potential, 3D-TDSE short range potential. The laser intensity 
is in unit of Jo with J 0 = 1 x 10 14 W/cm 2 , and the ionization 
potential is I p = 0.5 atomic units for all cases. 

tion of the terahertz field strength given by Ethz (t) oc 
X^E 2u (t) E* (t) E* (t) cos (j) and the OPD is close to 0 
or 7 r, in agreement with the numerical results at the 
laser intensity of Iq = 10 14 W/cm 2 . When TI domi¬ 
nate (7 < 1 ), the soft-recollision between ionized elec¬ 
tron with atomic core’s Coulomb potential play the key 
role [29, 30]. Rescattering currents dominate and es¬ 
caped currents vanish, and the OPD shifts from 7r to 
0 . 67 T. As further increasing of the laser intensity, OTBI 
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takes over for laser intensity over I = The non¬ 

vanishing initial longitudinal momentum cause the en¬ 
hancing of the escaped part of the electron wave packet, 
and hence the OPD is close to 0.57T. In this case, 
the atomic potential might be neglected and hence the 
photon-current model holds predicting the terahertz field 
as E THz (t) oc d J (t) /dt oc f \E U (t)] E 2u> (t) sin <j>. 

We now turn to the dependence of optimal terahertz 
yields (OTY) on the laser intensity. Surprisingly, the 3D 
calculation and the photocurrent model give quite sim¬ 
ilar dependence. It indicates that the ionization rate is 
crucial in determining the yields, since the ADK ioniza¬ 
tion rate is derived for 3D Hydrogen atom. In contrast, 
the ID and 2D calculations over-estimate the ionization 
yield, because the wave packet is less spreading compar¬ 
ing the 3D case. It is also worth noted that the yields 
from the short range potential are much lower for two 
reasons. First the tunneling ionization rate at the same 
ionization energy is indeed smaller for short range po¬ 
tential comparing to the long range potential. Secondly, 
the laser-assisted soft collision is suppressed in the case 
of short range potential. 

Next, we examine how the optimal phase delay and 
yields depend on the intensity ratio, a = I 2u /I u = 
E% u / E't , between the second harmonic and the funda¬ 
mental pulse. As shown in Fig. 3(a), the optimal phase 
delay remains nearly constant for a < 0.1 for given fun¬ 
damental laser intensity (/^ = 1.5 x 10 14 W/cm 2 ) of the 
fundamental pulse. But as the intensity ratio increases 
further, the OPD quickly drops to 0.57T, because the sec¬ 


ond harmonic pulse begins to participate in the ionization 
process which eventually take over the ionization from 
the fundamental pulse. It is shown in Fig. 3(b) that the 
OTY’s are scaled as power of the intensity ratio, which is 
consistent with the experiment [43]. The dependence on 
the intensity ratio of both the OPD and OTY are con¬ 
sistent with laser-assisted soft-collision model presented 
in [30]. For convenience of direct comparison to experi¬ 
ments, we present the intensity-ratio dependence by fix¬ 
ing the total power of the fundamental and the second 
harmonic pulse and the similar conclusion can be drawn 
as well. The ionization probability is also shown in Fig. 
3(b) which is less dependent on the intensity ratio, sug¬ 
gesting the dynamics following ionization is crucial in 
terahertz wave generation. 



FIG. 3. (color online) The scaling of OPD (a) and OTY 
(b) of THz generation with the intensity ratio of the second 
harmonic pulse to the fundamental laser pulse. 


To further clarify the mechanism of THz wave genera¬ 
tion by two-color laser pulses, we examine the wavelength 
scaling of THz generation. As the wavelength increases, 
the OPDs obtained from TDSE shift from 0.87T to 0.57T 
shown in Fig. 4(a), while the PC model predicts the 
OPD at 0.57T independent of wavelength. The variation 
of OPDs with wavelength might be rationalized as fol¬ 
lows. Because the wave packet is more diffused during 
each cycle for the longer wavelength, the soft-recollision 
with atomic nuclear plays less rule resulted in the OPD 
at 0.57T as predicted by PC model, which neglects the 
Coulomb potential completely. 

On the other hand, according to the PC model, the 
OTY’s are scaled quadratically with wavelength, while 
the TDSE simulation shows deviation from this depen¬ 
dence. As shown in Fig. 4(b), the OTYs from the TDSE 
calculation are scaled as A 1 ' 45 , at the laser intensity of 
Iu = 2/ 0 . While for I u = 1.5Jo or J 0 the deviation is 
even bigger. One of the reason for the deviation is the 
diffusion of the electron wave packet. However, the ex¬ 
perimental data from Clerici et al [44] is more close to 
the predication by PC model. This discrepancy needs 
further investigation in the future. 
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FIG. 4. (color online) Wavelength scaling of (a) the optimum 
phase delay and (b) THz yields. The experimental data is 
taken from [44]. Here, the two-color intensity ratio a = 1% 
and laser pulse has duration of FWHM 50fs. 


B. Terahertz wave generation from aligned 
molecules 


The previous investigation of THz generation from 
atoms in two-color laser fields demonstrates the impor¬ 
tance of Coulomb potential on the laser-driven electron 
dynamics. In addition it confirms the sensitivity of THz 
yields on the detuning of ionization and electron dynam¬ 
ics. For molecules, it is therefore expected that THz 
generation might exhibit dependence on both molecular 
alignment and molecular potential. 

To simulate molecular terahertz generation, we numer¬ 
ically solve the two-dimensional TDSE for a model H^ 
with the soft-core potential between the electron and the 
two nuclei, 


I/(r) 


1 

^/(r — R/2 ) 2 + a 2 


1 

•y/(r + R/2 ) 2 + a 2 ' 


where r = (z,x) and R = (R cos 8, R sin 8) are the elec¬ 
tronic and molecular coordinates, R is the internuclear 
distance. Two-color pulse is used. The fundamental and 
the second harmonic pluses are chose parallel in polar¬ 
ization along z axis. 8 is the alignment angle between 
the molecular axis and the laser polarization z , a is the 
soft-core parameter. Here, the motion of the nuclei are 
frozen and we set R = 2au, a = lau, and the ionization 
energy I p is equal to 0.934au. 

The emitted THz radiation are calculated from the low 
frequency part of the Fourier transformation of the in¬ 
duced dipole moments both parallel and perpendicular 
to the laser polarization. The total terahertz yield is 
the sum of yield of parallel G || and perpendicular Gj_ 
direction. As shown in figure 5(a), when the laser in¬ 
tensity is in the tunneling regime (7 > 1 , and / > 
2.1 x 10 14 W/cm 2 ), the OTY of terahertz field parallel 
to the laser polarization is two orders larger than that 
from the perpendicular direction for molecules aligned at 
45° with respect to the laser polarization. At 0° and 90° 
alignment, the perpendicular components vanish. There¬ 
fore the total THz yields are dominated by the emission 
of THz wave with field component parallel to the laser 
polarization and in the follows we make no distinguishing 


between the total yield and the yield parallel to the laser 
polarization. Figure 5(b) shows the OPDs of THz for 
molecules at different alignment. The OPD varies from 
7 r in the MPI range to 0.67T in the TI or OTBI range as 
the laser intensity increases, similar to the case of atoms. 




FIG. 5. (color online) Comparison of (a) the OTY and (b) the 
OPD of THz generation from molecules at different alignment 
and from atoms. 


For given two-color laser intensity, e.g. I u = 4 x 
10 14 W/cm 2 ,/ 2 oj = 0.5%/^,, we examine in details how 
the OTYs vary with alignment angle. Figure 6 (a) shows 
that the total OTYs take maximum at 8 = 0° and 
minimum at 8 = 90° respectivley, which closely assem¬ 
bles the angular dependence of the ionization probability. 
However, the perpendicular component of THz wave, al¬ 
though much weaker, exhibits very different alignment 
dependence, reaching maximum at 8 = 45° in contrast 
to the parallel component. 

We now investigate the alignment dependence of the 
OPDs. As shown in Fig. 6 (b), the OPDs of the total 
THz wave yields (mainly the parallel component) varies 
with the alignment angle 8. The maximum variation 
Ac/> m = (j) m (8) — (j) m (0) is found less than 0.17T. Al¬ 
though the variation is small, but it indeed indicates that 
the molecular potential plays a role in the generation of 
THz waves. In fact, converting into time scale, the phase 
change of 0.17T corresponds to a time delay of 67 attosec- 
onds (as). For better illustration, the variations A<f> m 
are plotted for both the parallel and perpendicular com¬ 
ponents in Fig. 6 (c). For both components, it can be 
seen that the OPD increases as the the alignment angle 
increases and reaches maximum at 90° alignment. The 
extra phase delay for molecules perpendicular to the laser 
polarization might be attributed to the larger recollision 
cross section in this geometry which leads to larger dis¬ 
tortion of the electron trajectories and hence more pro¬ 
nounced focusing effects during THz wave emission. 

In order to provide more insight into the generation of 
terahertz waves from molecules, we present the results of 
high harmonic generation (HHG) as well. In Fig. 7(b), 
the optimal yields of even harmonics from the two-color 
laser pulse are shown in contour plot against the align¬ 
ment angles. For comparison, high harmonic yields from 
the fundamental pulse alone are presented in Fig. 7(a). 
It can be seen that the yields of even harmonics gener¬ 
ated by the two color laser field peak at both alignment 
of 0° and 90°, which assembles the behavior of odd har- 
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FIG. 6. (Color online) Alignment-resolved OPD and OTY for 
THz generation from Hj". (a) Normalized OTY parallel (black 
solid) and perpendencular (red dash) to laser polarization, 
ionization population (blue dash-dot) and the 30th HHG (thin 
blue solid), (b) OPD. (c) OPD difference at various alignment 
angles to that of 9 = 0°. The laser pulse parameters are 
uj = 0.057au, Iui — 4/q, hu, = 0.5%/„, FWHM=15fs. 



FIG. 7. (Color online) Harmonic generation from H^". (a) 
Odd harmonic yields in one-color laser, (b) Even harmonic 
yields from the two-color laser pulse, (c) OPD difference of 
even harmonics relative to 9 = 0°. (d), (e), (f), the OPDs and 
optimum yields of harmonics 30th, 40th and 50th as functions 
of on alignment angles 6, respectively. 


monies generated by the one-color laser pulse. However 
this is in very contrast to the alignment dependence of 
terahertz radiation yields which peak at zero alignment 
(see Fig. 6(a) ). It is worth to note that the minimum 
of harmonic yields (shown in red lines in Fig. 7(a) and 
(b)) are in accordance to the predication (shown in black 
lines) that the destructive two-center interference occurs 
for 


i?cos 9 = (2m + 1) A/2, m = 0,1,..., (7) 

where A is the electron wavelength, and electron energy 
is determined by E/. = Noj [45]. 

By varying the phase-delay between the two pulses, 
both even harmonic and terahertz wave yields modulate. 
In Fig. 7(d-f), we show the OPD and OTY for 30th, 
40th and 50the harmonics as functions of alignment re¬ 
spectively. It can be seen wherever the OTY takes the 
minimum, the corresponding OPD makes a phase jump. 
In contrast, the OPD of THz wave show a slow variation 
with the molecular alignment (see Fig. 6(b)). 

For alignment of 0°, 45° and 90°, the OPD’s as func¬ 
tions of the harmonic order are shown in Fig. 8(a). It can 
be seen in all three cases that the OPD is nearly linear 
with respect to the harmonic order, reflecting the chirp 



Harmonic order 3 Alignment angle (0) Alignment angle (9) 


FIG. 8. (Color online) (a) The OPD of HHG at alignment 
angles of 0°, 45° and 90°. The black, red, and green dashed 
lines denote the OPD of terahertz waves at the corresponding 
angles. The inset shows the difference of OPD at 45° and 90° 
with respect to the alignment of 0° angle. The alignment- 
dependent (b) OPDs and (c) optimum yields of THz gener¬ 
ation and the yields of the 38th harmonics in the two-color 
laser pulse. 


of the recombine electron wave packet during HHG pro¬ 
cesses [46]. Detailed examination shows that the OPD 
for alignment 45° and 90° differs slightly from that of 0° 
alignment. 

To make more close contact of THz and HHG, the 
OPD and OTY of 38th harmonic yields are shown in com¬ 
parison with terahertz generation in Fig. 8(b) and (c). 
It clearly demonstrates that THz yields can be used to 
gauge the harmonic yield to gain the two-center interfer¬ 
ence information of HHG. The phase-delay dependence 
of THz generation can be used to probe the modulation 
of harmonics as well. 


IV. CONCLUSIONS 


In conclusion, we have investigated the atomic and 
molecular terahertz generation with intense two-coloar 
laser pulses. Based on the TDSE simulation, we show 
that THz generations is dominated by the competition 
between laser field and long-range coulomb field, which 
can affect the ionization mechanism and wave-packets 
dynamics. As the laser intensity increases, the OPD va- 
rations demonstrate the dominated mechanisms of THz 
generation are changing from four-wave-mixing, rescat¬ 
tering currents of soft-recollision between ionized elec¬ 
tron with atomic core, to photocurrent model without 
Coulomb potential. When we focus on THz genera¬ 
tion dependence on both molecular alignment angle and 
molecular potential, the total OTY is closely similar to 
the angular dependence of the ionization probability, and 
the OPD variation of several tens of attoseconds indicates 
that the molecular potential plays a role in the generation 
of THz waves. Finally we show that both the alignment- 
dependence of the OTY and OPD of THz generation, 
can be used to contrast that of high harmonic generation, 
which might be used to explore further the corresponding 
attosecond electron dynamics experimentally. 
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